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Abstract 

We work out the method for evaluating the QCD coupling constant 
at finite temperature (T) by making use of the finite T renormalization 
group equation up to the one-loop order on the basis of the background 
field method with the imaginary time formalism. The background field 
method, which maintains the explicit gauge invariance, provides noto- 
rious simplifications since one has to calculate only the renormalization 
constant of the background field gluon propagator. The results for the 
evolution of the QCD coupling constant at finite T reproduce partially 
the ones obtained in the literature. We discuss, in particular, the origin 
of the discrepancies between different calculations, such as the choice of 
gauge, the break-down of Lorentz invariance, imaginary versus real time 
formalism and the applicability of the Ward identities at finite T. 



*To appear in the Festschrift dedicated to Andrzej Bialas in honour of his 60th birthday, Acta 
Physica Polonica B27 (1996). 



1 Introduction 



One of the application of field theory at finite T [1-4] is to find the behaviour of 
coupling constant as a function of energy, temperature, and the chemical potential 
using the renormalization group (RG) equation. The knowledge of coupling constant 
at finite T environment then can be used, for instance, in perturbative calculations 
for quark-gluon plasma created in ion-ion collisions at high energies, in the evaluation 
of the grand unification scale in a cosmological context and the shift of the energy 
levels in a hydrogen plasma. 

A great deal of attention has been paid to the study of the behaviour of the QCD 
coupling contant at finite T [5-19]. The resulting formula for the temperature and 
scale dependent part of the coupling constant has proved to be very sensitive to the 
prescription chosen. In this paper we wish to reexamine this and clarify the origin 
of the discrepancies between different calculations in the literature. We employ the 
background field method (hereafter abbreviated as BFM) [20,21], which is based on 
a manifestly gauge invariant generating functional. The BFM provides notorious 
simplifications since we have to calculate only the renormalization constant of the 
background field (BF) gluon propagator. We discuss this method in section 2.1. 
Since the coupling constant depends not only on energy (//) but also on T, we derive 
a pair of RG equations - one for fi and one for T. This is done in section 2.2 where 
we also derive its solution. 

In order to get the renormalization constant of the BF gluon propagator we have to 
calculate the polarization tensor. At finite T due to the lack of the Lorentz invariance 
the structure of the polarization tensor is not equivalent to the one at T = 0. In 
section 3.1 we give a prescription how to define the renormalization constant at finite 
T. In section 3.2 we calculate the polarization tensor in the one-loop approximation. 
In section 3.3 calculating the vacuum part of the polarization tensor we reproduce 
the well known formula for T = [22,23], which reads as 

bV)]" 1 = [^oMrMl + 2g 2 (fi )K ln(fi/fi )}, (1-1) 

where 

K = (UN/3 - 2n F /3)/(47r) 2 . (1.2) 

Here N is for SU(N), rip is the number of flavours, fi is the energy scale and /x is 
the reference scale. 

We calculate the matter part of the polarization tensor in section 4.1 and, in sec- 
tion 4.2, we derive the formulas for T and fi dependent part of the coupling constant. 

In section 5.1 we review the results of the previous works and in section 5.2 we 
compare the asymptotic expansion formulas at a = ///T C 1 derived in different 
schemes. We discuss, in particular, the origin of the discrepancies between different 
calculations and give a few comments. 
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2 Formal methods 



2.1 Background field method 

Let us start with the generating functional for QCD: 

Z(J) = J DAD^D^Dr]Df]exp{i J d A x[C Q c D (A^ , f ) 
-(G a ) 2 /(20 + r a A^ a + fjdG a /du a V }}, (2.1) 

where A is the gluon field, \I/ and ^ are the fermion fields, r] and 77 are the ghost 
fields, (G a ) 2 is a gauge fixing term with £ being a gauge parameter, a; a is a S77(3) 
group parameter, and J Ma 7L M(l is a source term. In S77 (3) the gauge fields transform 
as: 

A^U(A^-id^/g)U\ (2.2) 
*->C/"*,t->W t , (2.3) 

V^Uri,fj^f}tf, (2-4) 
where [/ (x) = exp[iw a (a;)T a ] is a unitary transformation with T a being a generator for 
577(3). By using transformations (2.2)-(2.4) one can show that the gauge invariance 
of Z(J) is lost in commonly used gauges such as d^A^ = 0. Thus we see that that the 
explicit gauge invariance, which is present at the classical level in gauge field theories, 
is lost at the quantum level. 

The advantage of the BFM [20,21] is that it can maintain the explicit gauge 
invariance. For this purpose we divide the gluon field A^ into a sum of a classical BF 
Bfj, and a quantum field as 

A* = B » + ( 2 - 5 ) 

and choose for the gauge fixing term G a the BF gauge condition 

G a = (D,(B)) ab Ql (2.6) 
where is the covariant derivative: 

(D li ) ab = d li 6 ab + gf abc B°. (2.7) 
Then the generating functional reads 

Z(J,B) = J DQD^D^DrjDf] J exp{* J d A x[C QCD {B + Q, , f ) 
-(G°) 2 /(20 + J^A m + fjdG a /du a rj\}. (2.8) 
This functional can be shown to be gauge invariant by using the transformations: 

B^U{B^-idJg)U\ (2.9) 
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Q M -> t/Q^, (2.10) 

Dfi(B) — > UD^{B)U\ (2.11) 

J„(£) - UJ,{B)U\ (2.12) 

in addition to Eqs. (2.3) and (2.4). Actually we are only changing the integration 
variables in Eq. (2.8). 

In the BFM the quantum gauge fields and the ghost field need not be renormalized 
since they appear only inside loops. Thus only renormalizations 

9o = Z g g R , (2.13) 

B = {Z B f/ 2 B R) (2.14) 

£o = Z£ R , (2.15) 

are needed. 

The explicit gauge invariance of Z(J, B) implies that the perturbation series is 
gauge invariant in every order in g R and that Z B and Z g are related with each other. 
In order that the field tensor F£ v in Cqcd' 

FZ V = {Z B f' 2 [d,{Bl) R - d v {Bl) R + 9 R Z g (Z B )^r bc (B^) R (B^ R } (2.16) 
be gauge covariant one has to have 

Z g = (^)" 1/2 , (2.17) 

which enables us to calculate the evolution of coupling constant g R in a simple way. 
We assume this relation to be valid also at finite T. 



2.2 Coupled RG equations 

Equations (2.13) and (2.17) lead us to: 

g = g R (Z B y 1 / 2 . (2.18) 

Using the dimensional regularization [24,23] let us perform our calculations in (4 — 2e) 
dimensions. We notice from the QCD Lagrangian that the dimension of a gluon 
(quark) field is /i 1-e (/i 3 / 2-£ ,resp.). We redefine the bare coupling g so that it becomes 
dimensionless and rewrite Eq.(2.18) as 

go = g R (Z B )- 1/2 fi £ , (2.19) 

where g does not depend on T and /i. Then a pair of RG equations results upon 
taking the derivative of Eq. (2.19) with respect to T and \i 

Td/dT[g R {Z B )- 1 ' 2 ^] = 0, (2.20) 
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^d/d^[g R {Z B )- 1 / 2 ^} = 0. (2.21) 

A pair of these equations should determine the behaviour of with respect to the 
changes in T and \x. Generally Zg has the form [22,23] 

oo i 

Zb = 1 + E9r E4 5) /^ + / (!) (^)], (2.22) 

i=l j=l 

where A^/e^ are the divergent contributions independent on T and //, while /^(//, T) 
are convergent functions depending on T and /x which vanish in the low-T limit. 

Writing hereafter g = g^, for simplicity, one can derive the following RG equation 
from Eqs.(2.21) and (2.22) by taking the limit e -> 0: 

fidg/df, = g 3 {-AM + W2)d/diM[fW(iM,T)]} + 0(g 5 ). (2.23) 

This equation is reduced to the ordinary one-loop RG equation in the limit T — > 
since the function f^(fi,T = 0) vanishes. 

Similarly one can derive the RG equation for T from Eqs. (2.20) and (2.22) as: 

Tdg/dT = ( 9 y2)Td/dTifV(f,,T)] + 0(g 5 ). (2.24) 

Equations (2.23) and (2.24) constitute the coupled RG equations. One can easily 
obtain the solution to these coupled RG equations by integrating Eq. (2.23) from fj, 
to /j, and Eq.(2.24) from T to T. The desired solution reads as 

[g 2 (»,T)]-i = [g 2 (» ,T )}-i + 2A« In W/, ) 

-[/ (1) (^r)-/ (1) (Mo,T )]. (2.25) 

This equation describes the evolution of the QCD coupling constant as a function of 
T and \x [5-10]. (/io,T ) denotes the reference point. 



3 Polarization tensor 

3.1 Structure of the polarization tensor 

At T = environment the polarization tensor is Lorentz invariant and can be ex- 
pressed as [22] 

n^ = n(*V-W- (3.i) 

At T = it satisfies the transversality condition (current conservation): 

= o. (3.2) 
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At finite T, in the presence of matter, the Lorentz invariance is lost and the polariza- 
tion tensor can only be 0(3) rotational invariant [27,17]. Then it can generally depend 
only on 4 independent quantities, which we can choose, for example, Ilocfcjlloi and 
the two scalars Hl and Ut appearing in 

Uij = Yl T {5ij - kikj/k 2 ) + ULkikj/k 2 . (3.3) 

At finite T whether the transversality condition is satisfied or not depends on the 
gauge used. In the Coulomb gauge (9*^ = 0), for example, it is not transversal but 
is transversal at the one-loop level in the temporal axial gauge (A Q = 0), and in every 
order of the perturbative calculations in the BF gauge [17]. 
/,From Eq. (3.2) we have 

k U o = kiU i0 , (3.4) 

and from k^kyli^ = 0, we have 

^0-^00 = kikjUij. (3-5) 

Using Eqs. (3.3) and (3.5) we obtain 

n L = k 2 U 00 /k 2 . (3.6) 

For the coefficient of the transversal part of the polarization tensor one can derive 
from Eqs. (3.3) and (3.6) an expression 

n T = [n MM -n 00 (i + ^/fc 2 )]/2. (3.7) 

Thus we have derived in Eqs. (3.3)-(3.7) the general form of the T dependent 0(3) 
symmetric polarization tensor in the BFM. 

The polarization tensor can be splitted into a sum of a T independent (vacuum) 
part and a T and \x dependent (matter) part [4] : 

Since the polarization tensor at T = is related to the gluon propagator, D^ u , as 

= D R \ V - D~; u = UD^ U (3.9) 

and the relation between the bare and the renormalized propagators (see Eq. (2.14)) 
is 

IV = Z B l D RtlU) (3.10) 

one is led to 

Z B = l-U. (3.11) 

Thus the renormalization constant Z B can be obtained from the BF gluon self-energy 
tensor II. 
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In investigating the behaviour of the coupling constant at finite T one has to 
include the T dependent parts of II in Zb via Eq. (3.11). In determining Zb we 
encounter another type of ambiguity which is caused by the lack of the Lorentz 
invariance. In order to define II at finite T we generalize Eq.(3.11) as: 

%b — 1 — n|mc — n| matt , (3-12) 

where we have either 

n\ matt = n 00 /fc 2 , (3.13) 

or 

Tllmatt = n T /fc 2 . (3.14) 

n o and Ut are not connected with each other and in general there is no a priori way 
to decide which one is more natural. 



3.2 The polarization tensor at the one-loop level 

Our working formulas, which allow to analyze the evolution of QCD coupling constant 
at finite T, are: Eq.(2.22), Eq.(2.25) and Eqs.(3.12)-(3.14). Accordingly we have to 
evaluate the self-energy diagrams in Fig.l to get Z B up to the one- loop order. The 
Feynman rules for the interaction vertices are the same as in the T = case, and 
therefore identical to those given in [21]. Evaluating the one-loop diagrams (la)-(ld) 
in the Feynman gauge, we find for the boson contributions in the polarization tensor 

Urbason = W 2 N J d i p/(27r) 4 [Ag llu k 2 + 2{k il p v + k v p„) 

+Ap n p v - SkuK - 2(k + p) 2 g^]/{(k + p) 2 p 2 }. (3.15) 

For the fermion loop, neglecting the quark masses compared to /i and T, we have 
from the diagram (le): 

fermion = -4:ig 2 T F n F J d^pj (2ir) 4 {k^p u + k v p a 
+2puPu ~ 9uA(kp)+p 2 ]}/{(k + p) 2 p 2 }, (3.16) 

where T F = 1/2 for SU(3). 

Our polarization tensors satisfy the transversality condition 

ku^-/j,u\boson k fjTl ^j,] j ermion 0. (3.17) 
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3.3 Vacuum part of the polarization tensor 

Using a standard technique and introducing the Feynman parametrization in Eq. 
(3.15) we get the formula for the boson contributions: 

n Mi ,| 6oson = ig 2 N J d 4 p/(2n) 4 £ dx[Ag uu k 2 + 2(1 - 2x){k ti p v + k u p„) 
+^PnP» + (-3 -Ax + 4x 2 )k li k u ]/\p 2 + k 2 x(l - x)] 2 . (3.18) 

Similarly we get the formula for the fermion contributions from Eq. (3.16): 

Kfiulfermion = -2ig 2 n F J d 4 p/ (2n) 4 dx[2(-x + x 2 )k^k u + (1 - 2x) 
x(k u p u + k u p u ) - g uu {kp) + 2p ii p v + g uu k 2 x]/[p 2 + k 2 x{l - x)] 2 . (3.19) 

Notice that all of the integrals are ultraviolet divergent and thus have to be regu- 
larized. For this purpose we employ the dimensional regularization [24,23], which 
preserves gauge symmetries explicitly. The integrals for the vacuum parts become 
Euclidean if we change ip — > p 4 and thus can be easily evaluated. We obtain the 
following results: 

Urbason = -llg 2 N/(3a)(g uu k 2 - k„k v )/e + 0(1), (3.20) 

and 

IV I fermion 

2g 2 n F /(3a)(g uv k 2 - k n k u )/e + 0(1), (3.21) 

where 

a = (An) 2 . (3.22) 

Combining Eqs. (3.20) and (3.21) and taking into account Eqs.(2.22), (3.1) and (3.11), 
we reproduce the standard formula for defined in Eq. (2.25), i.e., Eq.(1.2) for 
T = 0. We end this subsection by emphasizing the simplicity of the calculation by 
the BFM in contrast to the conventional methods, e.g., in the covariant gauge [7-9]. 



4 Temperature dependent parts of the polariza- 
tion tensor 

4.1 Calculation of Iloo and Ilfifi 

To specify the subtraction point in calculating n o and we employ the static 
limit of zero external energy, which is commonly used in the literature [5]. In this 
prescription the momentum k is specified to be space-like k = (0, k) with k 2 = 
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— fi 2 . Such a choice enables us to determine the static properties. Making use of the 
imaginary time formalism [25,26] from Eqs. (3.15) and (3.16) we obtain for Hoo 

IIoo = g 2 T 2 (N/Q + n F /12)-2g 2 Nfi 2 [4F m (a)-F B2 (a)} (4.1) 
- 2g 2 n F/ 2 2 [F F o(a) - F F2 (a)}, (4.2) 

with a = (3fJL. Here the boson and fermion functions F in (/3fj,)(i =B,F) are defined as: 

F in (f3pi) = (1/a) / dxx n N l (^x/2)L, (4.3) 
Jo 

with 

L = In |(1 + x)/(l -x)\. (4.4) 
Similarly we derive an expression for H m : 

n w = g 2 T 2 (N/3 + n F /Q) - 2g 2 fi 2 [llNF B0 (a) + 2n F F F0 (a)}. (4.5) 



4.2 Temperature dependent part of the coupling constant 

As pointed out in section 3.1 we encounter an ambiguity in determining the renor- 
malization constant Z B as a direct consequence of the lack of Lorentz invariance. 
Here we write two formulas, one derived from n o and the other one from Hp. In the 
prescription with k = (0, k), one has from Eq.(3.7): 

n T = (n w - n 00 )/2. (4.6) 

Thus from Eq. (3.13) with n o we have 

-/ (1) (//, T) = (N/6 + n F /\2)/a 2 - 2N[AF m (a) - F B2 (a)} 
-2n F [F F0 (a) - F F2 (a)\, (4.7) 

and the other one from Eq.(3.14) with 

-/ (1) (/i, T) = [N/12 + n F /2A)/a 2 - N[7F B0 (a) + F B2 (a)} 
-n F [F F0 (a) + F F2 (a)}. (4.8) 
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5 Discussion 



5.1 Results of previous works 

Gendenshtein [5] calculated the QCD coupling constant at finite T in the one-loop 
approximation by using the RG equation, the dimensional regularization, and the 
covariant gauge with a space-like normalization momentum p M = (0,/x) and obtained 

-fV^,T) = (N/3 + n F /6)/a 2 . (5.1) 

Kajantie et al. [6] studied the gauge field part of QCD with N colors. They used 
the A = gauge and defined two renormalization schemes by writing 

D'l = D^(Z A - G/p 2 ) — (F — G)P^ (5.2) 

where D^ u is the gluon propagator and F, G and P^ v come from the polarization 
tensor 

IV = FP^ V + GPl„. (5.3) 

In the "magnetic prescription" they fixed the propagator at the point p M = (0, /i) and 
had for the T depending function the expansions (without quarks) 

- T) = 5JV/(16a), for a < 1, (5.4) 

- f (1) (fi, T) = iV[17/(90a 2 ) + 83a/(6300a 4 )], for a > 1. (5.5) 
In the so-called "electric prescription" they derived the coupling constant from F as 

- f^\fi,T) = N[l/(3a 2 ) - l/(4a) - 22(ln a)/ (3a)], for a < 1 (5.6) 

- / (1) (/i,T) = -iV[l/(18a 2 ) - lla/(900a 4 )], for a > 1. (5.7) 

Notice here the change of sign in the high-T behaviour. They concluded that the 
"magnetic prescription" is more natural than the "electric" one because the former 
uses the physical part of the gluon propagator. 

Nakaggawa et al.[7] used the real-time formalism and studied the scale-parameter 
ratios A(a)/A derived from different vertices in 4-flavour QCD. They used the co- 
variant gauge and found that the ratios derived from the three-gluon vertex and the 
gluon-ghost vertex show just the opposite behaviour than the one derived from the 
gluon-quark vertex. 

Fujimoto and Yamada [8] used the real-time formalism and derived the T depend- 
ing coupling constant from the gauge-independent Wilson loop. At a 1 , it reads 
as (without quark contributions) 

- f (1 \fi,T) = C]l/(3a 2 ) - l/(4a) - 22(lna)/(3a)]. (5.8) 
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The same authors have discussed the finite T RG equations [9] in the one-loop ap- 
proximation, using the real-time formalism, the covariant gauge and the dimensional 
regularization. Their results are summarized as follows. From the gluon propagator 
and three-gluon vertex: 

-/«(/!, T) = (C + T / )/(4a 2 ) + C[-23/3F B0 (a) - 3F B2 (a) - U/3G m (a) 
+32G B2 (a)} + 2T f [F F0 (a) - 3F F2 (a) + 152/9G F0 (a)], (5.9) 

where T f = 1/2, and C = N for SU(N). 

From the fermion propagator and fermion-gluon vertex: 

-fWfaT) = (C + 3C / + 2T / )/(12a 2 )-2(C / + llC)F B0 (a)/3 (5.10) 
- 16(67/ + 10C)G B0 (a)/9 - 32CG B2 (a) + 2(-5C/6 + C f /3 +(5/1)1) 
x F F0 (a) - 4(8(7/ + 17C)G F0 (a)/9 + lQCG F2 (a), (5.12) 

where C f = (N 2 - 1)/(2N) for SU(N). 

iFrom the ghost propagator and ghost-gluon vertex: 

-f (1 \fi, T) — (C + T/)/(12a 2 ) - C[7F m (a) + F B2 (a) 

+2G m {a)} - 2T f [F F0 (a) + F F2 (a)]. (5.13) 

The boson and fermion functions G in (a)(i =B,F) are defined as 

POO rl 

G in (a) = (2/ a) / dx / dyx^N^x^) /[x 2 (y 2 + 3) - 1]. (5.14) 

JO JO 

The functions G in (a)(i =B,F) do not appear in our results in Eqs. (4.6) and (4.7), 
since they come from the trigluon renormalization only. 

Stephens et al. [19] performed a BF one-loop calculation of gauge invariant beta 
functions at finite T, using the retarded/advanced formalism developped by Aurenche 
and Becherrawy [28] and derived: 

-fW(n,T) = (N/l2 + n F /24)/a 2 - N[21/AF m (a) 

+F B2 (a) + 7/2G B1 (a)} - n F [F F0 (a) + F F2 (a)}. (5.15) 

The numerical coefficients of the leading terms of the high-T expansion and the 
fermion parts are in complete agreement with our result in Eq.(4.7). In the boson 
parts, however, there are some numerical discrepancies with our result. 



5.2 Comparison of asymptotic expansions 

In order to compare our results with those mentioned in section 5.1 we derive the 
asymptotic expansions for —f^(fx,T). Equations (5.9)-(5.11) at a = fi/T <C 1 in the 
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high-T regime read as follows: 



-fW(^T) = nC/(9V3a 2 ) - 23C/(48a) + [(-22 

+7r/3v / 3)C + (18 - 387r/9v / 3)T / ](lna)/(3a), (5.16) 

-/ (1) (/i,T) = [(l-7r/v / 3)C + 3C , / ]/(12a 2 )-(llC , + C / )/(24a) 

-[(23 + 82^/9^)^ + (4 - 8n/9VS)C f + 47)] (In a) /(3a), (5.17) 

-/«(//, T) = -7C/(16a) - [C(22 + tt/VS) -87)] (In a)/(3a). (5.18) 
Next we derive the asymptotic expansion for Eqs. (4.6) at a <C 1. 

- (/i, T) = (JV/3 + n F /6)/a 2 - N/(2a) 

-[7V(224 - 71/3) - 4n/(£ 2 - 2/3)]/(3a), (5.19) 

where 

4 = ln(a/47r) + 7, (5.20) 

£ 2 = ln(a/vr) + 7. (5.21) 

The T 2 term coincides with the one of Gendenshtein [5], i.e.,Eq. (5.1). The gluon 
part in this asymptotic formula is consistent with the result of Elze et al. [17] derived 
in the BFM: 

- / (1) (/i,T) = A^[l/(3a 2 ) - l/(2a) -22(lna)/(3a) + ...], for a < 1. (5.22) 

This coincides with the result of Nadkharni [13], who has the terms up to 0(1/ a). 
Next from Eq. (4.7) we find the asymptotic formula: 

- / (1) (/i, T) = -7N/(lQa) - [7V(224 - 127/6) - 4n F (£ 2 - 5/12)]/(3a). (5.23) 

The behaviour of 11^ in the infrared region is known to have a form [17] 

limn r (0, k)//2 2 = g 2 [-cN/a + O(lna)]. (5.24) 

The factor c has been calculated in different gauges. In the covariant £-gauge its value 
is [27] 

c=(9 + 2£ + £ 2 )M (5.25) 
whereas in the temporal axial gauge it is [6] 

c = 5/16. (5.26) 

Our formula (5.21) at small a behaves as Eq. (5.22) with 

c = -7/16, for bosons, (5.27) 
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and 

c = 0, for fermions. (5.28) 

These results are consistent with the results of Refs. [9] [see Eq. (5.17)] and [13]. The 
relation (5.26) was also noticed by Elze et al. [17]. We note here that our c is negative 
and hence no spurious pole appears in the transversal propagator 

D Tlj (0, k) = -{5ij - hkj/k 2 )/^ 2 - IT T (0, fe)], (5.29) 

in contrast to the covariant and the temporal axial gauge cases. 

It is a known fact that the T 2 terms are gauge independent and the gauge param- 
eter dependence starts at ~ T order [14]. Thus the fact that the T 2 terms in Eqs. 
(5.14) and (5.15) are different from others, e.g. Eq. (5.1), should not originate from 
the gauge choices. 



5.3 Concluding remarks 

We end the paper by giving a few comments: 

1) In sections 5.1 and 5.2 we have seen that the T dependence of the QCD coupling 
constant is very sensitive to the prescription chosen. This is not 

a trivial issue, because all the results obtained hitherto also heavily depend on 
the vertex chosen (i.e., the trigluon, the ghost-gluon, or the quark-gluon vertex) to 
satisfy the renormalization condition of the QCD coupling constant. Furthermore in 
some gauges, e.g., in the Coulomb gauge, the transversality condition (3.2) on the 
polarization tensor does not hold and hence the structure of the polarization tensor 
in such gauges is different from the one which satisfies the condition. 

One of the reasons why we encounter different results in the literature is that 
the broken Lorentz invariance has not been treated as we have done in Eqs. (3.13) 
and (3.14). For example in deriving Eqs. (5.9)-(5.11) the authors of Ref.[9] have 
extracted the gluon and the vertex renormalization constants in front of Lorentz 
invariant structures not paying attention to the break-down of Lorentz invariance. 

Another possible source for discrepancies is that the Ward identities are used at 
finite T for different renormalization constants. The derivation of the Ward identity 
is based on the gauge invariance and also on the Lorentz invariance at T = 0. Thus 
it is not a surprise that results from different gauges or even (within the same gauge) 
from different vertices are totally different. 

The correspondence between the imaginary and real time formalisms has been 
investigated in detail [7,15-16]. The choice of the formalism can also be the source of 
the discrepancies under study was pointed out in these references. 

To clarify the issue of choosing a suitable renormalization prescription, one would 
need to compute the two-loop contributions to the coupling constant at finite T. It 
was shown, in particular, in the massive O(N) scalar model that the one-loop result is 
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drastically changed by two-loop contributions at high T and in zero momentum limit 
[29] . Recently Elmfors and Kobes have also argued that the inclusion of higher- loop 
effects is necessary to consistently obtain the leading term in the calculation of the 
thermal beta function in a hot Yang-Mills gas [30]. 

2) As we have seen the function / W (//, T) shows a power-like T-dependence instead 
of a logarithmic fall-off as a function of ji. Thus we have 

-\f {1) (n,T)-f m (lM>,T )] 

= £ c n [(T/f,r-(T /f, o r] + .... (5.30) 

n=l,2 

However, if a relation T = (const) \x holds, where (const) ~ 1/3 [6], then one would 
have a logarithmic T-dependence like in T = environment. Such a relation holds if 
one is dealing with the thermal equilibrium. 

We have derived Eq. (2.25) using two essentially different RG equations (2.23) and 
(2.24). Suppose we had used only one RG equation for T, i.e., Eq.(2.24), then we 
would have 

-[f {1 H»,T)-fW(» ,T)} 

= £ c n [(T/ri n -(T/n>) n ] + ..., (5-31) 

n=l,2 

instead of Eq.(5.28). In this case we would have a power-like dependence even with 
the relation T = (const) /i. 

We would like also to mention that all the machinary for the evolution of the 
running coupling constant at finite T can be analogously applied for the case of a 
quantum field theory (such as QCD) at finite energy as formulated in [31]. 

In conclusion, we note in accordance with the previous observations that there is 
in fact no unique way to define a T depending QCD coupling constant and the issue 
of finding its reasonable prescription is left as a subject of further investigation. 
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Figure Caption 

Fig.l. Diagrams for the one- loop calculation of the BF renormalization factor Zb'- 
a) Gluon loop; b) Ghost loop; c) Gluon loop from the 4-gluon vertex; d) Ghost loop 
from the 2-gluon 2-ghost vertex; e) Fermion loop. Wavy lines are quantum gauge 
field propagators. Wavy lines terminating in a "B" represent external gauge particles. 
Solid lines are fermion propagators and dashed lines represent ghost propagators. 
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